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ABSTRACT 

We report a serendipitous discovery of six very low-luminosity active galactic nuclei (AGNs) only by optical 
variability in one-month baseline. The detected flux variability is ~ 1-5% of the total luminosity of host galax- 
ies. Careful subtraction of host galaxy components in nuclear regions indicates that the fractional variability 
AF/F of the nuclei is of order unity. At least one of them is showing a compelling flaring activity within just 
a few days, which appears to be quite different from previously known AGN variability. We obtained spec- 
troscopic data for the one showing the largest flare and confirmed that it is in fact an AGN at z = 0.33 with an 
estimated black hole mass of ~ 10 8 M Q . As a possible interpretation, we suggest that these activities are com- 
ing from the region around the black hole event horizon, which is physically similar to the recently discovered 
near-infrared flares of our Galactic nucleus. It is indicated that our Galaxy is not special, and that surprisingly 
rapid flaring activity in optical/near-infrared bands may be commonly hidden in nuclei of apparently normal 
galaxies with low Eddington ratios, in contrast to the variability of well-studied luminous AGNs or quasars. 
Subject headings: black hole physics — galaxies: active 



1. INTRODUCTION 

Active galactic nuclei (AGNs) are generally variable, giving 
important information on the activity of central super massive 
black holes (SMBHs). In optical bands, most AGNs show 
significant variability on time scales longer than months, but 
day-scale or shorter variability is generally small (fractional 
amplitude of < 10%) and rare (Ulrich, Maraschi, & Urry 
1997; Webb & Malkan 2000; Hawkins 2002). Therefore the 
recent discovery (Genzel et al. 2003; Ghez et al. 2004) of 
near-infrared (NIR) flares within one hour from Sgr A*, an 
extremely low-luminosity AGN of our Galaxy with a SMBH 
mass of ~ 3 x 10 6 M Q , is surprising. Our Galaxy may be a 
peculiar AGN, or instead, such variability may be common 
for very low-luminosity AGNs like Sgr A*, but it could not 
have been examined because of the difficulty in observation 
for extragalactic sources. 

In a deep search for faint transient objects using the Sub- 
aru Telescope which has the largest field of view (30' x 24') 
among existing 8m class telescopes, we serendipitously found 
several AGNs only by variability on one month baseline, 
which is about 1-5% of the total luminosity of host galax- 
ies. There have been several attempts to find AGNs by their 
variability (Hawkins 1983; Bershady, Trevese, & Rron 1998; 
Trevese et al. 1989; Sarajedini et al. 2003), but most of such 
surveys looked for bright AGNs compared with host galaxy 
luminosity on time baselines of many (five to ten) years. 
Therefore our discovery is unique in terms of the time scale 
and nuclear/host luminosity ratios, providing us with a new 
method to find low-luminosity AGNs. Here we report this 
discovery and also suggest a possible (though not exclusive) 
interpretation. 
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2. THE OBSERVATION AND VARIABLE OBJECT SEARCH 

A fixed field centered on the galaxy cluster Abell 2152 was 
consistently monitored in V and / bands with similar exposure 
times in this observation, with the primary scientific purpose 
of searching microlensing events (Totani 2003). A reference 
frame was observed using a night on May 5, 2003, and the 
main-run observations during four consecutive nights were 
performed about one month after the reference observation 
(June 1-4). The weather condition was good and seeing was 
mostly stable; the effective seeing FWHM on the stacked im- 
ages over each night is 0.76, 0.68, 0.78, 0.56, and 0.62 arcsec 
for the reference frame and the four nights of the main run, re- 
spectively. After stacking these frames for each night, we sub- 
tracted the reference frame from the main-run frames by using 
the Alard & Lupton (1998) algorithm. We then searched for 
any brightened objects in subtracted frames (and in their in- 
verted frames for faded objects) by a standard software, SEx- 
tractor (Bertin & Arnouts 1996), requiring S/N > 10 of the 
SExtractor best magnitude both in V and / bands (about 26.1 
and 24.4 mag, respectively) in either of the four nights. The 
CCD counts were saturated in the centers of the brightest stars 
and galaxies, and these regions were removed from the anal- 
ysis because the subtraction does not work correctly. 

Using this procedure we detected more than 20 objects, 
and most of them are likely to be supernovae that are associ- 
ated with faint host galaxies but are offset from their centers. 
However, if we choose variable objects associated with bright 
galaxies (/ < 20), we found 6 variable objects, all of which are 
located at the very centers of well-resolved host galaxies (Fig. 

indicating that they are very likely to be AGNs. Assuming 
that supernovae trace stellar light, the chance probability of 
finding a supernova in the center of a host galaxy within the 
angular resolution is about 3%, and hence the probability that 
all these six objects are supernovae by chance is ~ 2 x 10~ 10 . 
The objects 1—4 brightened in one month, but the other two 
faded. 

To estimate the statistical significance of these detections, 
we measured an aperture flux (1.2" = 6 pixels diameter) of 
centers of all ~ 1,000 galaxies with / < 20 on the subtracted 
images. The distribution is fit by a Gaussian, and the S/N ra- 
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FIG. 1 . — V and / color images of the six variable AGNs found in the survey. The first column shows the reference images taken about one month prior to the 
main observing run, and the reference-subtracted images of the consecutive four nights in the main run are shown from the second to fifth columns. The size of 
each panel is 10" X 10", and the pixel scale is 0".2/pixel. The upper four objects (1—4) brightened in a month, while the lower two faded. The subtracted images 
of the lower two are inverted for presentation. The V and / surface brightness profiles of host galaxies are also shown in the right panels; photon counts are in 
arbitrary unit and two-dimensional surface brightness data are projected onto the primary axis of elliptical fitting. 



TABLE 1 

Properties of Variable AGNs and Thier Hosts 



V band / band 



ID 






p 


A/w 


AF/F [%\ 




^V.tot 




p 


A/v, c 


AF/F L%J 


I'd L*J 


1 


32.10 


3.18 


-/0. 15/0.21 


0.34 (9.2<t) 


-/227/162 


84 


87.98 


12.79 


0.70/1.79/2.20 


1.50(13(7) 


214/84/68 


11 


2 


22.87 


4.59 


0.44/0.81/0.96 


1.41 (34cr) 


320/174/147 


0.0 


72.15 


20.47 


3.58/5.47/6.38 


1.71 (14<t) 


48/31/27 


7.9 


3 


9.92 


3.15 


0.20/0.58/1.02 


0.81 (21<x) 


405/140/79 


67 


36.45 


13.88 


2.84/4.42/6.04 


0.68 (6.2a) 


24/15/11 


4.2 


4 


15.83 


3.47 


0.05/0.45/0.53 


0.44(12(7) 


880/98/83 


0.0013 


71.52 


17.95 


0.90/3.26/4.15 


1.26 (llcr) 


140/39/30 


6.7 


5 


9.20 


2.39 


0.25/0.54/0.60 


0.13 (3.8(7) 


52/24/22 


46 


41.51 


14.89 


3.43/4.92/5.00 


1.48 (13.8ct) 


43/30/30 


29 


6 


19.83 


4.96 


0.27/0.89/1.56 


0.36 (llcr) 


133/40/23 


23 


69.00 


19.95 


3.27/5.73/7.74 


0.34 (3.6ct) 


10/5.9/4.4 


46 



NOTE. — The host galaxy total flux (F wAol ), nuclear fluxes within 1.2" diameter aperture in the reference frame (F UmC ), estimates of AGN components after star-light subtraction 
(7W-*)> variability flux (AF UX -, one month baseline, maximum among the four nights), and the fractional variability (AF/F = AF W , c / F u . c -*) are shown. All fluxes are in /jJy. The 
normal photometric error within the aperture is ~ 0.01 and 0.04 /jJy for V and / bands, respectively, while the errors used to calculate S/N of AF VmC (shown in parentheses) include 
those associated with the image subtraction procedure. The three numbers separated by slashes for F UjC -* and AF/F are those assuming three different profiles of host galaxy light: 
the r lj/4 law, the exponential profile, and constant surface brightness within the aperture, respectively. For V-band data of the object 1, star- subtracted nuclear flux is negative when the 
r^ 4 -law profile is used, and hence it is not shown. Significance of day-scale variability is indicated by P$, which is the chance probability of getting the data from a constant source 
during the four nights. The objects 1—4 brightened from May to June, but the other two faded. 
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FIG. 2. — (Left panels) Variability flux (in /dy) of the six AGNs during 
the four nights of the main observing run, relative to the reference frame 
taken one month earlier. The filled and open circles are for V and / bands, 
respectively. (Right panels) Spectral index a (/'„ oc u a ) of variability flux 
between V and / bands. The error bars in the left panel and thick error bars in 
the right panel are statistical 1 a errors, while the thin error bars in the right 
panel are the total including systematics (see text). The upper four objects 
brightened in a month, but the lower two faded and their flux is measured in 
inverted subtracted images. 



tios of the maximum variability flux within the same aperture 
among the four nights are estimated by using this empirical 
noise level. The results are shown in Table 1. All objects 
are statistically significant using a criterion of S/N > 10 in 
at least one band. Three of the six objects show this month- 
like timescale variability at such high levels of significance in 
both bands. In principle, we cannot exclude a possibility that 
these are false detections created by some systematic failure 
of image subtraction. However, failures of subtraction, which 
often happen in saturated regions, can easily be recognized by 
complicated and clearly artificial image profiles. The satura- 
tion effect is sensitive to seeing, and if these detections were 
created by saturation, the variability flux of all objects should 
show a similar dependence on the observed date, which is not 
actually observed. Therefore we consider that these detec- 
tions of one-month scale variability are reliable. As shown 
below, the spectroscopy of the object 2 shows some AGN fea- 
tures, which provides a further support for the significance of 
the variability detection. 

3. FLUX VARIABILITY ESTIMATES 

The fluxes of the variable objects measured on the sub- 
tracted images during the four days are shown as light curves 
in Fig. |2] The host galaxy luminosities are shown in Tabled 
and the variability flux is typically about 1-5% of these. We 



estimated detection efficiency and photometric errors by plac- 
ing artificial point sources on the surface brightness peaks of 
the ^1,000 galaxies with / < 20. Then we repeated the same 
procedure of image subtraction and source detection. We 
found that the detection efficiency is about 50% for /„ =0. 176 
and 0.515 /iJy, or 25.8 and 24.1 mag for V and / bands, re- 
spectively. 

We made two different error estimates: one is "statistical", 
which is the variance from the mean of the observed flux of 
an artificial point source during the four nights, and the other 
is "total" including systematics, which is the variance from 
the real flux of an artificial point source. The significance of 
variability within the four days should be estimated by the 
statistical errors. The chance probability of getting the data 
from a constant source during the four days is shown as P4 
for each object in Tabled The evidence for variability within 
the four days is compelling for the object 2 and strong for the 
object 4. The spectral indices between V and / bands are also 
shown in Fig. |2] with the total error estimates. The evidence 
of spectral variability within days is weak for all objects. The 
values of the indices show a large scatter from object to object, 
but strong conclusions cannot be derived because of the large 
systematic errors. 

The AGN component in the pre-subtraction images must 
be estimated to know the fractional variability ( AF /F) of the 
nuclei. There is no apparent central excess in surface bright- 
ness profiles of host galaxies (Fig. [0, indicating that stellar 
light dominates AGN luminosity. We estimate the host galaxy 
component in the nuclear regions of the reference frame as 
follows. We set an aperture diameter of I", 2 as the nuclear 
region. The surface brightness profiles outside this aperture 
are elliptically fitted by the r'/ 4 -law as well as the exponen- 
tial profiles after convolved with the seeing. Then we subtract 
them from the observed total nuclear fluxes, to get an esti- 
mate of the AGN components in the reference frame (Table 
0. The use of the exponential profile is conservative because 
it gives lower host galaxy contribution and hence lower frac- 
tional variability estimates, compared with the r 1//4 -law that 
is often used for the central stellar light profile. To be even 
more conservative, we also estimate the host galaxy contribu- 
tion by a constant surface brightness within the aperture, with 
a value estimated by photon counts along the aperture annu- 
lus. The dominant error is that of variability flux AF, which 
can be inferred from S/N ratios in the table. These results 
indicate fractional variabilities (AF/F) of order unity for the 
six objects. 

4. SPECTROSCOPY AND REDSHIFT ESTIMATES/DETERMINATION 

Although there are two galaxy clusters at z = 0.04 and 0. 13 
in the observed field (Totani 2003), the locations of the six 
AGNs are not concentrated to the cluster centers. The mor- 
phology of these host galaxies is apparently elliptical or early- 
type, with magnitude of / ~ 18.5-19.5. The field galaxy 
counts at / ~ 19 dominate the expected number of the clus- 
ter galaxies. According to a galaxy evolution model that is 
in good agreement with various data (Totani & Yoshii 2000), 
the peak of the redshift distribution of early type galaxies at 
this magnitude is z ~ 0.3-0.4, and observed V -I colors (~ 
1.6-2.2) are consistent with expected values for passively- 
evolving early-type galaxies in this redshift range. Therefore, 
it is likely that these galaxies are field early-type galaxies at 
z ~ 0.3-0.4, i.e., out of the two clusters. 

We obtained an optical spectroscopic data by Sub- 
aru/FOCAS for the object 2 (Fig. [3}, and confirmed that 
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FIG. 3. — The optical spectrum of the object 2. The redshift is z = 0.33. 



the redshift is, as expected, z = 0.33 and the host galaxy is 
a typical giant galaxy with Mb ~ -20. Several emission lines 
are detected in the nuclear region, whose FWHMs (~ 300 
km/s) are typical of narrow line regions. The line ratios are 
typical of the low ionization nuclear emission regions (LIN- 
ERs), confirming that this object is a low-luminosity AGN. 
In the following, we assume z = 0.4 for the other five AGNs 
for order-of-magnitude discussions. Using the empirical rela- 
tion between bulge luminosity and black hole mass (McLure 
& Dunlop 2001), the SMBH mass of these AGNs should be 
about ~ 10 8 M Q , and the nuclear luminosity becomes ~ 10~ 5 - 
10~ 4 in units of the Eddington ratio, i.e., (yL v ) / 'Lead, where 
Lndd is the Eddington luminosity. 

5. DISCUSSION 

For the first time, except for the recently discovered flares 
in Sgr A*, fractional AGN variability of ~ 100% within a few 
days in optical/infrared bands is detected in one, and prob- 
ably six, such low-luminosity (low Eddington ratio) AGNs. 
The Sgr A* flares have been interpreted as violent nonthermal 
phenomena at the innermost region of a radiatively inefficient 
accretion flow (RIAF) around the SMBH (Yuan, Quataert, & 
Narayan 2004). Such behavior cannot be seen for luminous 
AGNs whose time variability has been studied much better; 
the only exception is the rare population of blazars whose jets 
are closely directed to the observer (Ulrich et al. 1997). Al- 
though yet other interpretations might be possible, we here 
examine these two interpretations (RIAFs and blazars). 

The continuum-to-line luminosity ratio is a useful diagnos- 
tic to discriminate between the two interpretations. Since the 
continuum of blazars is strongly beamed by relativistic jets, 
unbeamed emission lines are generally weak and even dif- 
ficult to see. On the other hand, if we are observing less 
beamed continuum emission from accretion disks, we expect 
a continuum-to-line ratio that is typical of normal (non-blazar) 
AGNs. The Ha emission (2.6 x 10 40 erg/s) is clearly detected 



for the object 2 and it argues against the blazar interpreta- 
tion. By using the empirical relations between line and con- 
tinuum luminosities (Ho & Peng 2001; Wang, Staubert, & Ho 
2002), the Ha luminosity translates into the expected opti- 
cal continuum luminosity of vf v ~ 1.6x 10 41 and 7.8 x 10 43 
erg/s for normal AGNs (Seyfert Is) and blazars, respectively. 
This should be compared with the variability flux of the ob- 
ject 2: A(uf v ) ~ (0.9-3.2) xlO 42 erg/s in the V band. The 
observed variability flux is intermediate between the two in- 
terpretations, but considering that the object 2 is likely in a 
flare phase, the mean flux appears to be closer to the RIAF 
interpretation. 

The number density of these six AGNs is about ~ 
10~ 4 Mpc~ 3 , using the comoving volume out to z = 0.4. If the 
AGN variability is from disk emission, this number can be 
compared with published AGN luminosity functions in op- 
tical (Ho 2004) or X-ray (Ueda et al. 2003) bands at a nu- 
clear activity level of vL v ~ 10 42 erg/s, which is found to be 
4>(> 10 42 erg/s) ~ 10~ 35 Mpc" 3 . These numbers are not very 
different, indicating that the six AGNs found in this work con- 
stitute a considerable fraction of low-luminosity AGNs. On 
the other hand, if variability is from blazar activity, the real 
number density should be much larger than ~ 10 -4 Mpc" 3 , be- 
cause of the collimation of blazar emission. It is believed that 
blazars are beamed by a factor of Ail /(4ir) ~ T" 2 /4 < 1 /400, 
where r ~ 10 is the jet Lorentz factor (Salvati, Spada, & 
Pacini 1998), and hence the real number density of AGNs like 
the objects found in this work should be > 10 Mpc -3 . How- 
ever, the isotropic Ha line luminosity indicates that the nu- 
clear activity must be larger than ~ 10 41 erg/s, and the number 
density of such AGNs should not be larger than ~ 10~ 3 Mpc" 3 
from the AGN luminosity functions. This provides another 
argument against the blazar interpretation. 

It is interesting to note that, if the variability reflects the 
scale of emission region around SMBHs, the variability time 
scale of ~ 1 hour for the 3 x 10 6 M© black hole of Sgr A* cor- 
responds to ~ 1 day for ~ 10 8 M Q black holes. Although the 
Eddington ratio of our AGNs (~ 10~ 5 -10~ 4 ) is much larger 
than that of the Sgr A* (~ 10 -9 ), it is well below the bor- 
der (~ 10~ 2 ) separating the standard thin disk and RIAFs 
(Kato, Fukue, & Mineshige 1998). A RIAF model of NIR/X- 
ray flares of Sgr A* predicts that flare activity should dis- 
appear with increasing Eddington ratio, but it remains up 
to L/L Edd ~ 10" 5 in optical/NIR bands (Yuan et al. 2004). 
Therefore it seems a reasonable interpretation that previously 
unknown violent activity in optical bands has been hidden 
around the SMBH event horizon of very low Eddington-ratio 
AGNs, which cannot be seen in luminous AGNs whose lumi- 
nosity is close to the Eddington limit. 
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